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1.  introduction 


In  this  paper,  we  model  the  time  evolution  of  a  point  on  the  centerline  of  a  dispersing  plume 
from  a  solid  fuel  rocket  booster  (in  particular,  a  Titan  IV  rocket)  in  the  stratosphere.  We  will  do  this 
at  two  altitudes — 20  km  and  30  km — ^to  look  for  the  transient  effect  on  the  stratospheric  ozone  layer. 
At  this  writing,  there  are  some  other  papers  in  the  literature  (Refs.  1, 2, 3, 4)  modeling  this 
phenomenon,  but  the  earlier  papers  do  not  include  the  effects  of  chlorine  oxide  dimer  reactions  on  the 
ozone,  and  the  Denison,  et  al.,  paper  is  for  a  much  smaller  rocket  and  uses  a  smaller  set  of  chemical 
reactions.  Our  reaction  set  is  based  on  a  survey  of  several  literature  models  of  stratospheric 
chemistry,  and  we  believe  that  the  gas  phase  reaction  set  is  complete  to  the  extent  of  present 
knowledge.  Our  starting  point  will  be  the  conditions  predicted  by  the  JANAF  plume  code  at  2 
seconds  after  passage  of  the  rocket  (see  Zittel,  Ref.  5).  At  this  time,  there  will  be  a  considerable 
amount  of  free  Cl  and  CI2  from  afterburning  of  the  plume  in  the  atmosphere.  Figure  1,  quoted  from 
Zittel,  shows  the  fraction  of  HCl  converted  to  active  form  a.s  a  ftmction  of  altitude.  This  figure  is 
based  on  the  chemistry  taking  place  0-2  seconds  after  passage  of  the  exit  plane.  The  plume 
temperatures  covered  by  the  Zittel  model  range  from  about  2,000  K  down  to  400  K.  We  will 
compute  the  subsequent  (2  seconds  to  10  days)  (400  K  to  ambient  temperature)  altitude  dependent 
reactions  of  chlorine,  nitrogen,  and  hydrogen  oxides  that  occur  in  the  reactively  mixing  plume.  We 
will  then  discuss  the  implications  of  the  model  results  and  give  some  conclusions  and 
recommendations  for  increasing  our  understanding  of  this  transient  phenomenon. 


2.  Background 


Hydrogen  chloride  (HCl)  is  a  major  combustion  product  of  solid  rocket  boosters  that  use 
ammonium  perchlorate  oxidizer.  As  this  gas  leaves  the  motor,  plume  models  run  at  Aerospace  and 
elsewhere  (Refs.  1,  5)  indicate  that  afterburning  takes  place,  converting  a  substantial  amount  of  the 
HCl  to  Cl  and  CI2  (21-5%,  depending  on  altitude).  This  occurs  through  reactions  such  as: 

OH  +  HCl  =  Cl  +  H2O  (1) 


H  +  HCl  =  Cl  +  H2 


(2) 


The  chlorine  so  released  will  immediately  react  with  ozone  to  begin  a  catalytic  destruction  cycle.  At 
high  altitudes  (30  km,  for  example),  the  sequence  will  be: 


Cl  +  O3  =  CIO  +  O2 


(3) 


O  +  CIO  =  Cl  +  O2 


(4) 
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At  low  altitudes  (e.g.,  20  km),  a  more  complex  ozone  destruction  cycle  may  take  place 


involving  the  chlorine  oxide  dimer  (C10)2  (Ref.  6): 

Cl  +  O3  =  CIO  +  O2  (3) 

CIO  +  CIO  =  (C10)2  (5) 

Cl  +  (C10)2  =  CI2  +  ClOO  (6) 

ClOO  +  M  =  Cl  +  O2  (7) 

CI2  +  hv  =  2C1  (8) 


This  cycle  is  similar  to  that  which  causes  the  Antarctic  ozone  hole,  except  that  the  chlorine 
oxide  dimer  forms  because  of  the  high  concentrations  of  chlorine  in  the  rocket  plume,  and  not  because 
of  low  temperatures.  Also,  as  we  discovered  in  ranning  this  model,  the  principal  route  for  generation 
of  atomic  chlorine  from  the  dimer  is  not  photolysis  of  the  dimer  (which  is  the  route  in  the  Antarctic) 
but  the  branched  chain  reaction  of  the  dimer  with  atomic  chlorine  as  shown  above.  Since  the  two 
reaction  cycles  shown  above  regenerate  atomic  chlorine,  these  cycles  will  continue  to  destroy  ozone 
until  the  Cl  and  CIO  return  to  HCl  by  reactions  with  hydrogen-containing  molecules  such  as  methane 
or  hydrogen: 


Cl  +  CH4  =  HCl  +  CH3  (9) 

In  the  normal  stratosphere,  the  Cl  and  CIO  are  typically  converted  to  HCl  by  this  process  in  a  few 
days,  with  the  HCl  returning  only  very  slowly  to  active  form.  In  a  rocket  plume,  however,  there  is 
sufficient  Cl  and  CIO  so  that  substantial  local  depletion  of  ozone  may  take  place  on  shorter  time 
scales.  There  is  some  early  observational  evidence  that  a  local  ozone  hole  may  be  created  (Ref.  7). 

There  is  a  very  slow  mixing  of  the  stratosphere  with  the  troposphere,  so  that  over  a  period  of 
several  years,  HCl  released  into  the  stratosphere  will  eventually  return  to  the  troposphere  and  rain  out, 
with  an  average  characteristic  time  of  a  few  years. 
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3.  Comparison  with  Chlorofluorocarbons 


During  the  time  it  takes  for  HCl  to  diffuse  back  to  the  troposphere,  it  will  spread  over  the 
globe,  and  one  must  consider  how  HCl  from  rockets  will  compare  with  HCl  from 
chlorofluorocarbons  (CFCs).  A  typical  CFC  has  a  very  long  lifetime  in  the  troposphere,  and  since  the 
troposphere  contains  about  95%  of  the  mass  of  the  atmosphere,  95%  of  a  typical  CFC  resides  passively 
in  the  troposphere.  At  any  given  time,  the  5%  present  in  the  stratosphere  is  slowly  photolyzing  into 
atomic  chlorine: 


CFC  +  UV  light  =  Cl  +  organic  radicals  (10) 

The  chlorine  atoms  produced  by  photolysis  will  enter  ozone  destruction  cycles  as  described 
above.  Because  CFCs  have  been  produced  for  many  years,  the  total  amount  in  the  atmosphere  is  on 
the  order  of  25  million  tons,  and  about  300,000  tons  of  atomic  chlorine  are  produced  each  year  from 
photolysis.  Because  of  the  long  lifetime  of  CFCs,  this  production  of  atomic  chlorine  is  expected  to 
continue  with  a  characteristic  lifetime  of  about  60  years  even  in  the  absence  of  further  manufacture  and 
release  of  these  molecules.  In  comparing  stratospheric  chlorine  release  from  rockets  with  that  from 
CFCs,  the  global  effect  is  thus  similar  to  a  hypothetical  CFC  with  a  lifetime  of  a  few  years  (Ref.  19). 


4.  Potential  Heterogeneous  Reactions 


Neither  the  global  nor  local  effects  described  above  consider  the  potential  additional  effects  of 
the  particulate  exhaust  of  the  solid  rocket  boosters.  At  this  time,  any  heterogeneous  chemistry  that 
takes  place  on  the  particles  (mostly  alumina  and  aluminum  oxychloride)  is  poorly  understood.  Some 
reactions  that  have  been  proposed  are: 

1 .  the  direct  destraction  of  ozone  by  surface  catalysis: 

2O3  =  /via  particles/ 30  2  (11) 

2.  catalysis  of  the  photolysis  of  CFCs: 

CCI3F  +  2H2O  +  hv  =  /via  particles/ CO  2  +  HF  +  3HC1  (12) 

3.  catalysis  of  reactions  of  “reservoir”  species  to  put  chlorine  back  into  more  reactive 
form,  e.g.: 


CINO3  +  HCl  =  /via  particles/ Cl  2  +  HNO3 


(13) 
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Reactions  such  as  this  are  analogous  to  the  reactions  that  take  place  on  polar  stratospheric  clouds  to 
produce  the  Antarctic  ozone  hole,  but  they  are  just  beginning  to  be  examined  in  the  laboratory  (see 
reference  K  at  the  end  of  the  kinetics  data  table).  The  potential  effect  of  such  chemistry  is  one  of  the 
things  we  will  examine. 


5.  Model  Description 


The  computer  model  is  based  on  the  SURFACE  CHEMKIN  code  developed  by  R.  J.  Kee  of 
Sandia,  Livermore  (Refs.  8, 9, 10).  This  package  contains  subroutines  for  calculating  thermodynamic 
properties,  equilibrium  constants,  rates,  and  unit  conversions,  for  parsing  character  data  such  as  species 
and  element  names,  and  reaction  mechanisms,  and  for  solving  the  set  of  differential  equations. 
SURFACE  CHEMKIN  is  a  set  of  software  tools  and  a  subroutine  library  designed  to  help  create  a 
program  to  solve  a  specific  chemical  kinetics  problem.  The  package  is  designed  to  solve  a  group  of 
differential  equations  subject  to  a  set  of  constraints;  various  properties  of  the  system  may  be  extracted 
along  the  way.  One  needs  to  write  the  program  which  describes  the  chemical  system  of  interest,  and 
calls  the  appropriate  subroutines  to  do  the  calculations.  Multiple  surface  and  bulk  phases  can  be 
included.  Stuface  rates  are  scaled  by  the  ratio  of  surface  area  to  gas  phase  volume.  The  package  is 
designed  to  model  chemical  vapor  deposition.  We  wrote  a  driver  program  incorporating  some  changes 
in  the  calculation  method,  added  some  new  subroutines,  and  modified  some  of  the  original  routines. 
The  resulting  program  is  well  adapted  to  atmospheric  chemistry.  Our  alterations  are  described  in  more 
detail  below. 

The  plume  model  is  adiabatic.  CHEMKIN  calculates  the  enthalpy  of  each  species  in  the 
model.  As  species  are  created  or  destroyed  through  reaction  they  add  or  subtract  energy  from  the 
plume;  mixing  with  the  air  also  effects  the  energy  balance.  This  net  change  in  energy  is  divided  by  the 
heat  capacity  of  the  plume  to  obtain  the  temperature  change  for  each  integration  step.  This  energy 
conservation  equation  is  included  as  a  governing  equation,  or  constraint,  for  SURFACE  CHEMKIN. 

The  governing  equation  for  temperature,  alluded  to  above,  is  based  on  conservation  of  energy, 
and  is  included  here  for  reference. 


N 

i=l 

ra«f  ^  ^  1 

dt  dt  ^  dt  ^ 

+  1 

i^Air  _ 

1  pAir 

1 

CpP  + 

) 

The  preceding  equation  describes  the  change  in  the  plume  temperature  with  time.  The  numerator  in  this 
equation  is  the  total  enthalpy  change  for  the  plume  in  time  interval  dt .  The  summation  is  over  species 

/ ,  where  N  is  the  total  number  of  species  considered;  here  N  is  25.  is  the  enthalpy  in  mass  units, 
m,-  is  the  molecular  weight  of  species  i,  and  is  the  number  of  moles  produced  by  gas,  g,  surface. 


s,  or  liquid  phase,  i,  reactions  per  unit  volume  (area  for  phase  =  ^);  the  product 


dnf 

dt 


represents 
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the  heat  gained  or  lost  to  the  plume  due  to  the  production  (or  consumption)  of  a  particular  chemical 
species  in  gas  phase  reactions  in  the  time  interval  dt .  The  rate  of  species  production  from  surface 

reactions,  ■  is  normalized  by  the  ratio  of  surface  area  to  gas  phase  volume.  is  the  ratio  of 
dt 

liquid  phase  volume  to  gas  phase  volume,  which  normalizes  the  rate  of  species  production  from  bulk 

phase  reactions,  3^.  Tis  the  temperature,  Cp  is  the  cross-sectional  area  of  the  plume,  and  is 
dt 

heat  capacity.  The  term  -  t)  — ^  represents  the  heat  required  to  bring  ambient  air,  which 

is  mixing  into  the  plume,  to  the  same  temperature  as  the  plume  gases.  The  superscript  Air  refers  to 
quantities  in  the  undisturbed  air  outside  the  plume.  The  denominator  in  the  above  equation  represents 
the  specific  heat  of  the  plume  and  is  made  up  of  the  specific  heat  of  the  gas  phase,  C|p,  and  condensed 

phase,  Cpm£,  portions,  where  p  is  the  gas  phase  density,  and  the  total  mass  in  the  liquid  phase  per  unit 
volume  of  gas  is  .  The  heat  capacity  of  the  surface  phase  is  negligible  and  is  ignored. 

The  remaining  constraints  are  based  on  conservation  of  mass  and  are  broken  down  into 
separate  equations  for  each  species  in  each  phase.  The  governing  equation  for  gas  phase  species  is: 


dYi 


dt 


+ 


Cpdt 


The  symbol  Y  is  used  for  the  gas  phase  mass  fraction.  The  first  term  in  this  equation  is  the  molar 
changes  due  to  reaction  converted  to  a  mass  fraction,  and  the  second  term  is  the  effect  of  dilution  on 
plume  species.  This  equation  has  a  similar  stmeture  to  the  numerator  of  the  first  governing  equation. 

ITie  equation  for  surface  species  in  terms  of  site  fraction,  X,,  is  similar.  There  is  also  a 
correction  for  coverage,  where  the  density  of  surface  sites  in  a  given  surface  phase  j  is  given  by  T j  and 
the  number  of  sites  occupied  by  a  given  species  is  Oj-  Large  molecules  could  occupy  several  sites,  but 
all  the  species  in  this  model  are  small  molecules  with  a,-  =  1  •  Only  one  surface  phase  (alumina/gas)  is 
used  in  our  model.  There  is  an  equation  like  the  following  one  for  each  surface  species. 


dXi 


dt 


+ 


V«J 


This  complete  set  of  constraints  is  used  in  solving  the  differential  equations  from  the  kinetics 
for  each  time  step  in  the  simulation. 


To  model  the  alumina  particles  we  use  a  fixed  condensed  phase  volume  and  surface  area. 
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6.  Input  Data  for  the  Computations 


The  model  is  formulated  as  a  two  bin  system,  consisting  of  a  plume  bin  that  slowly  mixes 
with  the  surrounding  stratosphere,  and  a  stratosphere  bin  that  is  assumed  to  be  infinite,  i.e., 
concentrations  of  species  in  the  overall  stratosphere  are  assumed  not  to  be  perturbed  by  the  plume. 
Interbin  transport  is  described  by  a  simple  dilution  function  which  gives  the  centerline  concentration  as 
a  polynomial  function  of  time.  The  dispersion  (dilution)  of  the  plume  appears  to  be  the  most  poorly 
understood  part  of  this  phenomenon  at  this  time.  We  studied  several  different  rates  of  dispersion  to  get 
an  idea  of  the  coupling  between  the  chemistry  and  the  dispersion,  but  most  of  the  data  presented  is 
based  on  a  dispersion  model  published  by  JPL  (Ref.  11).  This  dispersion  model  was  fit  to  a 
concentration  that  was  a  polynomial  function  of  time.  In  all  the  mns,  we  include  0.01%  helium  in  the 
initial  plume  as  a  “tracer”  to  show  the  rate  of  dilution.  (The  natural  abundance  of  helium  is  omitted 
from  the  stratospheric  bin  for  simplicity.) 

The  model  includes  34  chemical  species  and  some  100  gas  phase  chemical  and  photochemical 
reactions,  and  two  heterogeneous  reactions  shown  in  Table  1  (see  page  1 1).  All  the  important 
stratospheric  reactions  of  nitrogen,  hydrogen,  oxygen,  and  chlorine  containing  molecules  are  included. 
Heterogeneous  reactions,  including  direct  loss  of  ozone  on  alumina,  and  catalysis  of  the  reaction  of 
chlorine  nitrate  with  HCl,  are  included  as  described  in  the  results.  Details  of  the  parameters  for  this  ate 
given  in  Appendix  A.  The  starting  conditions  at  20  km  and  30  km  are  given  in  Figures  2  and  3, 
respectively.  At  20  km  altitude,  the  initial  plume  temperature  is  400  K  and  the  surrounding 
stratosphere  is  217  K.  At  30  km  altitude,  the  initial  plume  temperature  is  438  K  and  the  surrounding 
stratosphere  is  227  K.  The  starting  conditions  are  based  on  a  JANAF  plume  model  (Ref.  5)  for  a  Titan 
rv  sized  launch  vehicle  (Ref.  12).  The  hole  sizes  are  based  on  the  hydrogen  chloride  release  rate  of  2.0 
tons  per  kilometer  at  20  km  altitude  and  1 .8  tons  per  kilometer  at  30  km  altitude  (Ref.  6). 


7.  Description  of  Results 


The  model  results  are  shown  as  log  of  the  species  concentrations  versus  log  of  the  time,  and 
temperature  versus  log  of  the  time.  Figure  4  shows  the  results  at  20  km  and  Figure  5  the  results  at  30 
km.  These  two  figures  are  for  the  complete  model,  including  heterogeneous  chemistry, 
photochemistry,  and  chlorine  and  nitrogen  oxide  chemistry.  It  can  be  seen  that  the  concentration  of 
ozone  is  reduced  by  about  a  factor  of  100  for  about  lO"^  seconds  at  30  km  altitude,  and  for  about  2  x 
10  seconds  at  20  km.  In  each  figure,  the  final  ozone  concentration  is  the  ambient  concentration  and  the 
final  temperature  is  the  ambient  temperature. 

Several  features  common  to  both  curves  deserve  comment.  First,  for  the  initial  hundred 
seconds  less  than  1%  of  the  ozone  remains.  There  is  a  partial  resurgence  of  the  ozone  between  100  and 
1,000  seconds.  This  appears  to  be  related  to  the  fall  in  temperature  at  this  time  as  the  cooling  effect  of 
stratospheric  mixing  is  felt.  There  is  also  a  pronounced  peak  in  the  OCIO  and  in  the  chlorine  oxide 
dimer  concentration  after  10“*  seconds  at  both  altitudes.  This  is  due  to  the  recovery  in  the  ozone 
concentration,  which  rapidly  reacts  with  any  remaining  Cl.  The  dis^pearance  of  the  Cl  removes  a 
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major  sink  for  the  dimer  and  so  there  is  a  resurgence  in  the  concentration  of  the  dimer  between  10  and 
10^  seconds  before  final  dispersal  in  the  stratosphere. 

Figure  6  shows  the  results  at  20  km  altitude  with  the  photolysis  reactions  “turned  off.”  The 
primary  effect  here  is  that  there  are  a  few  Cl  atoms  from  the  afterburning  (see  Figure  1),  but  no  new 
formation  of  Cl  atom;  thus  Cl  2  accumulates.  There  is  no  catalytic  cycle  here;  each  Cl  atom  can 
destroy  one  ozone  molecule  at  most  before  being  ultimately  tied  up  as  (C10)2  or  as  CI2.  This  stops 
the  ozone  destruction  chemistry  after  100  seconds,  and  the  ozone  recovers  due  to  mixing  with  the 
ambient  stratosphere  beyond  this  time. 

Figure  7  shows  the  results  at  20  km  altitude  with  heterogeneous  chemistry  and  NO^  chemistry 
eliminated.  The  overall  results  are  similar  to  the  full  reaction  mechanism,  as  shown  in  Figure  4, 
indicating  the  dominance  of  gas  phase  chlorine  chemistry  in  the  plume.  This  is  not  surprising  given 
that  there  is  10^  times  as  much  chlorine  in  the  plume  as  in  the  ambient  stratosphere,  whereas  there  is 
only  10“*  times  as  much  NOj^.  The  surface  reactions  require  either  ozone,  which  is  already  depressed 
by  Ae  gas  phase  reactions,  or  CIONO2,  which  is  present  at  low  concentrations  in  the  ambient 
stratosphere  and  mixes  slowly  into  the  plume.  Surface  reactions  may  be  important  at  longer  times  or  at 
the  plume  edges  where  gas  phase  chlorine  concentrations  are  smaller. 

Figure  8  shows  the  results  at  20  km  altitude  with  no  initial  chlorine  in  the  plume.  The  only 
chlorine  species  appearing  in  the  calculations  come  from  the  ambient  stratosphere.  There  is  only  slight 
ozone  depletion  due  to  NO^  chemistry,  heterogeneous  chemistry,  and  thermal  decomposition.  A 
comparison  with  Figure  7  confirms  that  gas  phase  chlorine  chemistry  dominates  ozone  loss  in  the 
plume.  Interestingly,  a  comparison  with  Figure  6  shows  the  tremendous  damage  done  by  the  chlorine 
atoms  initially  present  in  the  plume  during  the  first  100  seconds. 

Figure  9  shows  the  results  at  20  km  altitude  with  no  initial  NO^  or  chlorine  in  the  plume.  This 
is  very  similar  to  Figure  8,  and  shows  that  the  small  amount  of  ozone  depletion  seen  is  due  to 
heterogeneous  chemistry  and  thermal  decomposition;  NO^  does  not  contribute  significantly  to  ozone 
depletion  in  the  plume. 

Figure  10  shows  the  results  at  20  km  altitude  with  all  of  the  Cl  initially  in  the  form  of  HCl, 
i.e.,  with  no  “afterburning”  of  HCl.  Interestingly,  there  is  still  some  ozone  depletion,  although  not 
nearly  as  much  as  in  the  afterburning  case.  Even  beyond  the  afterburning  region  the  OH  in  the  plume 
is  still  much  higher  than  in  the  ambient  stratosphere.  We  attribute  the  ozone  loss  in  this  region  to 
“cold”  afterburning  (i.e.,  below  400  K).  (Note  that  there  is  still  some  atomic  chlorine  seen  in  this 
figure  due  to  reaction  of  the  HCl  with  OH  in  the  plume.) 

Figure  1 1  shows  double  curves  with,  and  without,  heterogeneous  chemistry,  and  with 
“enhanced”  heterogeneous  chemistry  (see  Table  1).  It  can  be  seen  that  including  the  heterogeneous 
chemistry  can  widen  and  deepen  the  ozone  hole  somewhat,  but  this  requires  a  heterogeneous  reactivity 
about  100  times  higher  than  that  seen  in  our  laboratory  experiments. 


7 


Figure  12  shows  the  effect  of  plume  dispersion  rate  on  the  persistence  of  the  “hole.”  For  these 
runs,  a  simple  Gaussian  dispersion  in  which  the  cross-sectional  area  of  the  plume  increases  as  the 
square  of  time  is  used.  Thus,  the  concentration  of  helium  decreases  with  the  square  of  time  multiplied 
by  a  parameter.  For  our  main  runs,  a  polynomial  in  time  was  used  (see  above).  It  appears  that  the 
slower  the  dispersal,  the  longer  the  hole  persists.  There  is  very  little  experimental  data  on  plume 
dispersion  rates. 

Figure  13  shows  the  approximate  peak  hole  size  as  a  function  of  dispersion  rate,  as  described 
for  Fig.  12.  The  hole  size  in  this  case  is  based  on  the  total  chlorine  concentration  at  the  time  the  ozone 
recovers,  using  the  HCl  loadings  from  a  Titan  IV  rocket  (Ref.  12).  The  size  of  the  hole  depends  to 
some  degree  on  what  is  meant  by  ozone  recovery;  in  this  analysis,  ozone  is  considered  to  recover  when 
it  reaches  30%  of  the  ambient  concentration. 


8.  Conclusions 


The  results  of  the  detailed  model  calculations  generally  support  our  initial  back-of-the-envelope 
calculations  on  the  size  of  transient  ozone  holes  created  by  passage  of  a  solid  rocket  booster  engine 
(Ref.  6).  Since  the  plume  dissipation  rate  is  the  most  poorly  understood  part  of  this  phenomenon,  we 
treated  this  rate  as  a  parameter  to  be  varied  over  a  wide  range  of  values.  The  model  has  given  the 
following  information  to  date:  1)  The  peak  hole  size  increases  if  the  plume  dissipates  more  slowly.  The 
slowest  rates  of  dissipation,  which  appear  to  be  favored  by  the  JPL  model  (Ref.  11),  give  a  hole 
roughly  10  km  in  diameter  at  20  km  altitude,  and  28  km  in  diameter  at  30  km  altitude.  This  result  is 
somewhat  counter-intuitive,  and  appears  to  be  due  to  the  second  order  nature  of  the  chemistry  involved. 
Since  the  formation  of  chlorine  oxide  dimer  will  depend  on  the  square  of  the  active  chlorine 
concentration,  the  chain  length  for  ozone  destraction  apparently  also  depends  on  the  square  of  the 
concentration.  Thus,  a  rapid  dissipation  of  the  plume  will  result  in  shorter  chain  lengths,  and  less 
overall  ozone  destmction.  2)  Similarly,  the  slower  the  plume  dissipates,  the  longer  the  hole  persists, 
lasting  from  a  few  hours  to  a  day.  These  model  results  suggest  that  dilution  is  the  primary  factor 
determining  the  effect  of  chlorine  on  ozone  depletion.  3)  Another  interesting  result  is  the  importance  of 
including  chlorine  oxide  dimer,  CI2O2,  in  the  reaction  scheme  at  lower  altitudes  such  as  20  km.  This 
species,  which  is  of  key  importance  in  the  formation  of  the  Antarctic  ozone  hole,  was  omitted  from 
most  previous  models  of  rocket  plume  chemistry.  Inclusion  of  this  species  makes  a  significant 
difference  in  the  hole  size  at  20  km  altitude.  4)  In  addition,  the  model  shows  that  photolysis  of  this 
species — neither  the  cross  section  nor  the  products  are  well  established  experimentally — (Refs.  20, 21, 
22, 23)  is  not  critical  for  ozone  destraction  in  the  first  few  hours  (it  may  be  important  at  later  times  or 
near  the  edges  of  the  plume).  In  the  special  situation  found  in  the  plume  center,  it  is  reaction  of  the 
dimer  with  atomic  chlorine  that  sustains  the  ozone  destraction  cycle. 
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9.  Recommendations 


The  transient  perturbation  of  the  stratosphere  by  a  passing  rocket  provides  a  unique 
opportunity  to  better  our  understanding  of  stratospheric  chemistry,  and  to  test  our  knowledge  and 
theories.  This  preliminary  work  suggests  that  an  observable  local  ozone  depletion  should  take  place, 
and  studies  of  this  phenomenon  are  important  for  understanding  the  environmental  impact  of  launch 
vehicles  and  for  basic  knowledge  of  the  atmosphere.  We  make  the  following  recommendations  for 
further  work: 

1 .  Sensitivity  studies  should  be  done  to  determine  the  gas  phase  rate  coefficients  most 
critical  to  determining  the  size  and  depth  of  the  local  ozone  hole,  and  which  need 
further  laboratory  study.  These  reactions  may  not  be  the  same  as  those  determining 
the  steady  state  ozone  concentration  in  the  global  stratosphere. 

2.  Similarly  new  heterogeneous  reactions  should  be  tested  for  importance  in  this 
phenomenon.  Thus  far  only  a  handful  of  such  processes  have  been  considered. 

3.  New  modeling  studies  of  the  potential  global  impact  of  rocket  launches  need  to  be 
done  in  view  of  our  changing  understanding  of  global  ozone  depletion.  In  particular, 
we  refer  to  the  newly  discovered  importance  of  heterogeneous  reactions  in  the 
natural  sulfate  aerosol  layer  (the  Junge  layer),  and  the  importance  of  redistribution  of 
ozone  depleted  air  from  the  Antarctic  (Ref.  24).  These  effects  account  for  a 
significant  change  in  the  predictions  of  current  models,  and  may  change  the  potential 
impact  of  rocket  launches  on  the  ozone  layer  as  well. 

4.  Planning  in  situ  and  remote  observations  of  this  phenomenon  must  be  done  in  concert 
with  the  predictions  of  a  good  model.  Some  spectroscopic  studies  by  Beiting  (Ref. 
18)  have  already  shown  the  importance  of  basing  observational  capabilities  on 
realistic  concentration  predictions  for  the  various  species  and  their  interferences. 
Further  work  is  necessary  for  each  new  observational  technique  that  is  proposed. 

5.  The  chemistry  and  photochemistry  of  chlorine  oxide  dimer,  CI2O2,  is  of  central 
importance  to  this  phenomenon,  as  well  as  the  Antarctic  ozone  hole.  In  our  case,  the 
photolysis  is  important  at  the  edges  of  the  hole,  and  at  the  time  the  hole  begins  to 
disappear.  The  chemical  decomposition  of  the  dimer  is  important  at  all  times,  and  is 
not  well  documented. 

6.  The  rate  of  mixing  of  the  rocket  plume  with  the  surrounding  atmosphere  is  poorly 
understood,  yet  the  model  calculations  of  ozone  depletion  are  sensitive  to  this 
parameter.  Measurements  are  needed  to  determine  the  rate  of  mixing  and  dissipation 
of  the  rocket  plume. 

We  will  be  glad  to  collaborate  with  interested  parties  in  developing  work  plans  for  pursuing 
each  of  these  research  topics. 
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Table  1.  Kinetic  Data 


Elements 

H  O  N  A1  Cl  C  He 

Species 

02  N2  O  0(1d)  oh  nos  H02  H202  Cl  CIO  CHS 
H20  NO  N02  OS  HCl  ClOO  OCIO  HOCl  C12  C1202 
HNOS  N204  C10N02  H  N  N20  H2  HNO  CO  C02 
H0N0N205CH4He 

Reactions 

*-E3/RT  3-1-1  6 

k  =  Ai  e  k  units  =  cm  molecule  sec  ,  Ea  in  calories,  for  3-body  reactions,  k  in  cni  molecule 

2-1  1 

sec  .  Photochemical  reactions  are  in  sec  ,  at  20  km/SO  km,  as  quoted  in  Reference  B  or  calculated 
from  the  data  therein.  (Local  noon,  April  1, 40  N.)  Photochemical  reactions  must  be  put  into 
computer  as  irreversible,  or  the  reverse  reactions  will  be  unrealistically  speeded  up  by  the 
thermodynamic  program. 


Reactants 

Products 

A 

n 

Ea(cal) 

References 

03  +  C1 

= 

010  +  02 

2.9E-11 

0 

520 

B 

OH  +  C12 

= 

H001  +  01 

1.4E-12 

0 

1,790 

B 

0  +  HOC1 

= 

OH +  010 

1.0E-11 

0 

4,370 

B 

C1  +  02  +  M 

=> 

0100 +  M 

1.4E-29 

-1.5 

0 

B 

C1  +0100 

= 

012  +  02 

2.3E-10 

0 

0 

B 

01  +0100 

= 

010  +  010 

1.2E-11 

0 

0 

B 

01  +  H02 

= 

HOI  +  02 

1.8E-11 

0 

-340 

B 

01  +  H02 

= 

OH  +  01O 

4.1E-11 

0 

890 

B 

010  +  H02 

H001  +  02 

4.8E-13 

0 

-1,390 

B 

010  +  0 

= 

01  +02 

3.0E-11 

0 

-140 

B 

0  +  012 

= 

010  +  01 

4.2E-12 

0 

2,720 

G 

01  +  H202 

= 

HOI  +  H02 

1.1E-11 

0 

1,950 

A 

0  +  0100 

= 

010  +  02 

1.4E-10 

0 

0 

H 

OH +  0100 

= 

H001  +  02 

1.4E-10 

0 

0 

H 

H  +  0100 

= 

HOI  +  02 

1.4E-10 

0 

0 

H 

O  +  O  +  M 

02 +  M 

5.2e-35 

0 

-1,790 

0 

H02  +  NO 

= 

N02  +  OH 

3.7E-12 

0 

-477 

A 

0  +  N02 

NO +  02 

6,5E-12 

0 

-238 

A 

N02  +  H 

= 

NO  +  OH 

4.0E-10 

0 

675 

B 

NO  +  0 

= 

N  +  02 

3.5E-15 

1 

38,747 

D 

0  +  N2 

= 

NO  +  N 

1.3E-10 

0 

75,506 

E 

N20  +  M 

= 

N2  +  0  +  M 

1.2E-6 

-0.73 

62,789 

D 

0  +  N20 

N2  +  02 

1.7E-10 

0 

28,017 

D 

N  +  N02 

= 

N20  +  0 

3E-12 

0 

0 

B 

N20  +  H 

N2  +  0H 

1.6E-10 

0 

15,101 

D 

N02  +  H2 

= 

HONO  +  H 

3.98E-11 

0 

28,998 

1 

HNO  +  H 

= 

H2  +  N0 

3E-11 

0 

994 

D 

HNO  +  OH 

H20  +  NO 

8E-11 

0 

994 

D 

NO  +  M 

= 

N  +  O  +  M 

2.4E-9 

0 

148,429 

D 
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Reactants 

Products 

A 

n 

Ea(cal) 

References 

N2  +  M 

= 

N  +  N  +  M 

3.82E5 

-3.5 

225,008 

E 

N  +  OH 

=: 

NO+H 

3.8E-11 

0 

-169 

A 

CIO  +  NO 

= 

N02  +  C1 

6.4E-12 

0 

-576 

B 

HOC1  +  OH 

H20  +  C10 

3E-12 

0 

994 

B 

H2  +  OH 

= 

H20  +  H 

5.5E-12 

0 

3,974 

B 

OH  +  0 

zz 

H  +  02 

2.2E-11 

0 

-238 

B 

OH  +  H 

= 

H2  +  0 

8.1E-21 

2.8 

3,875 

C 

OH  +  OH 

=: 

H20  +  0 

4.2E-12 

0 

477 

B 

OH +  00 

= 

C02  +  H 

1.5E-13 

0 

0 

A 

03  +  0 

= 

02  +  02 

8E-12 

0 

4,093 

B 

02  +  hv 

=> 

0  +  0 

2.0E-13/3.0E-11 

0 

0 

B 

0  +  02  +  M 

03 +  M 

3.0E-28 

-2.3 

0 

B 

03  +  hv 

=> 

0  +  02 

6.0E-4/6.0E-4 

0 

0 

B 

03  +  hv 

=> 

0(1D)  +  02 

3.0E-5/1.5E-4 

0 

0 

B 

0(1D)  +  H20 

= 

20H 

2.2E-10 

0 

0 

A 

0(1D)  +  M 

= 

0  +  M 

2E-11 

0 

-213 

A 

03  +  H02 

= 

OH  +  02  +  02 

1.4E-14 

0 

994 

B 

OH +  03 

H02  +  02 

1.6E-12 

0 

1,867 

B 

OH  +  OH  +  M 

H202  +  M 

6.6E-29 

-0.8 

0 

B 

H  +  OH  +  M 

= 

H20+M 

6.11  E-26 

-2.0 

0 

C 

H  +  02  +  M 

=: 

H02  +  M 

5.2E-28 

-1.6 

0 

B 

H02  +  OH 

H20  +  02 

4.8E-11 

0 

-500 

B 

OH  +  H202 

s= 

H20  +  H02 

2.9E-12 

0 

318 

B 

H  +  03 

=r 

OH +  02 

1.4E-10 

0 

954 

B 

H  +  H02 

OH  +  OH 

2.8E-10 

0 

874 

C 

H  +  H02 

= 

H2  +  02 

1. IE-10 

0 

2,130 

C 

H02  +  H02 

=r 

H202  +  02 

2.2E-13 

0 

-1,190 

A 

0  +  H02 

=: 

OH +  02 

2.7E-11 

0 

-445 

A 

0  +  H202 

OH  +  H02 

1.4E-12 

0 

3,970 

A 

03  + NO 

= 

02  +  N02 

1.8E-12 

0 

2,720 

A 

N03  +  0 

= 

N02  +  02 

IE-11 

0 

0 

B 

N03  +  NO 

= 

N02  +  N02 

1.8E-11 

0 

-219 

A 

N03  +  OH 

= 

N02  +  H02 

2.3E-11 

0 

0 

A 

N03  +  H02 

= 

N02  +  OH  +  02 

4.3E-12 

0 

0 

A 

0C10  +  C1 

= 

C10  +  C10 

3.4E-11 

0 

-318 

B 

C12  +  hv 

=> 

Cl  +C1 

2.0E-3/2.0E-3 

0 

0 

B 

C10  +  hv 

=> 

Cl  +0 

1.0E-4/5.0E-4 

0 

0 

B 

C100  +  hv 

=> 

C10  +  0 

.01 /.01 

0 

0 

B 

0C10  +  hv 

=> 

C10  +  0 

.3/.3 

0 

0 

B 

H0C1  +  hv 

=> 

0H  +  C1 

4E-4/5E-4 

0 

0 

B 

C1  +  HOC1 

HC1  +C10 

3.0E-13 

0 

0 

H 

C10  +  C10  +  M 

C1202  +  M 

8.7E-23 

-3.9 

0 

B 

CIO +  03 

= 

C100  +  02 

IE-17 

0 

0 

B 

C10  +  03 

= 

0C10  +  02 

IE-12 

0 

7,948 

B 

OH  +  HC1 

= 

H20  +  C1 

2.4E-12 

0 

656 

A 

0H  +  0C10 

H0C1  +  02 

4.5E-13 

0 

-1,590 

A 

OH  +  C100 

H02  +  C10 

1E-12 

0 

0 

H 
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Reactants 

Products 

A 

n 

Ea(cal) 

References 

H202  +  C10 

= 

HOC1  +  H02 

1E-13 

0 

0 

H 

C1202  +  C1 

= 

C12  +  C100 

1E-10 

0 

0 

A 

C1202  +  hv 

=> 

C100  +  C1 

1,3E-3/1.7E-3 

0 

0 

B 

C10  +  N02  +  M 

= 

C10M02  +  M 

4.8E-23 

-3.4 

0 

B 

C10N02  +  C1 

= 

C12  +  N03 

6.8E-12 

0 

-318 

A 

C10N02  +  hv 

=> 

C1  +  NOS 

6.0E-5/8.0E-5 

0 

0 

B 

C10N02  +  0 

C10  +  N03 

3.0E-12 

0 

1,589 

A 

C10N02  +  0H 

HOC1  +  NOS 

1.2E-12 

0 

655 

A 

C1  +  CH4 

= 

CHS  +  HCl 

9.6E-12 

0 

2,682.5 

A 

C1  +H2 

= 

H  +  HC1 

3.7E-11 

0 

4,570.1 

A 

0  +  0C10 

= 

C10  +  02 

2.5E-12 

0 

1,888 

A 

OC10+NO 

= 

N02  +  C10 

2.5e-12 

0 

1,192 

A 

HC1  +0 

:= 

OH  +  C1 

1.0E-11 

0 

6,636 

A 

H  +  C12 

= 

HC1  +C1 

1.43E-10 

0 

1,172 

G 

2H  +  M 

H2  +  M 

1.5E-29 

-1.3 

0 

C 

2C1  +M 

= 

C12+M 

6.14E-34 

0 

-1,800 

G 

H  +  O-f  M 

= 

OH  +  M 

1.3E-29 

-1 

0 

C 

CO  +  0  +  M 

= 

C02  +  M 

1.7E-33 

0 

3,000 

C 

NO  +  0  +  M 

= 

N02  +  M 

9.2E-28 

-1.6 

0 

A 

OH  +  N02  +  M 

HN03  +  M 

2.6E-22 

-3.2 

0 

B 

H  +  NO  +  M 

HNO  +  M 

2.47E-28 

-1.3 

735 

D 

OH  +  NO  +  M 

= 

HONO  +  M 

6.5E-25 

-2.4 

0 

A 

N02  +  03 

= 

NOS  +  02 

1.2E-13 

0 

4,868 

A 

H202  +  hv 

=> 

OH  +  OH 

1.0E-5/1.5E-5 

0 

0 

B 

H02  +  hv 

=> 

H  +  02 

1.0E-6/3.0E-6 

0 

0 

B 

N02  +  hv 

=> 

NO  +  0 

1.0E-2/1.0E-2 

0 

0 

B 

NOS  +  hv 

=> 

N02  +  0 

0.3/0.3 

0 

0 

B 

NOS  +  hv 

=> 

NO +  02 

0.03/0.03 

0 

0 

A 

HN03  +  hv 

=> 

OH  +  N02 

1.0E-6/1.0E-5 

0 

0 

B 

Heterogeneous  Reactions 

203/surface/  =  3O2  Y  =  10"^  normal,  10'^  “enhanced”*  J 


CINO3  +  HCl  =  HNO3  +  CI2  Y=  0.02  normal,  0.2  “enhanced”*  K 


*  See  description  for  Figure  1 1  in  text. 
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Figure  1.  Far- field  Chemical  Identity  of  Chlorine  as  a  Function  of  Altitude 
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Figure  2.  Starting  Conditions — 20  km 
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Figure  3.  Starting  Conditions — 30  km 
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Figure  4.  Ozone  Depletion  at  20  km,  Full  Reaction  Mechanism 
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Figure  5.  Ozone  Depletion  at  30  km,  Full  Reaction  Mechanism 
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Figure  6.  Ozone  Depletion  at  20  km,  No  Photolysis 
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Figure  7.  Ozone  Depletion  at  20  km,  Gas  Phase  Chlorine  Reactions  Only 
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Figure  8.  Ozone  Depletion  at  20  km,  NO,  and  Surface  Only 
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Figure  9.  Ozone  Depletion  at  20  km,  Surface  Reactions  Only 
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Figure  10.  Ozone  Depletion  at  20  km,  All  Chlorine  as  HCl,  No  NO, 
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Figure  11.  Ozone  Depletion  at  20  km  at  Three  Sticking  Coefficient  for  Surface  Reactions 
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Figure  14.  Trimodal  Size 
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Appendix  A 


Estimate  of  the  Surface  Area  of  the  Particles 

The  heterogeneous  reactivity  of  the  alumina  particles  will  be  proportional  to  the  surface  area 
available  per  unit  volume  of  air  and  the  sticking  coefficient  for  the  reaction  in  question,  according  to 
the  standard  kinetic  theory  result: 


dn/dt  =  nvYA/4 


^  _ 

where  n  is  the  concentration  of  the  reacting  species  in  molecules/cm  ,  v  is  the  molecular  velocity,  y 
is  the  sticking  coefficient,  and  A  is  the  available  surface  area  per  unit  volume  in  cm^/cm^.  The  area, 
in  turn,  will  be  determined  by  the  size  distribution  of  the  alumina  particles  and  their  shape.  For  our 
purposes,  these  properties  may  be  summarized  as  an  available  surface  area  per  gram  of  the  alumina, 
which  can  then  be  put  into  our  model.  We  know  the  mass  of  alumina  released  per  unit  altitude  for  the 
various  launch  systems  (Ref.  12),  and  this  can  be  combined  with  a  dispersion  model  to  give  an 
absolute  value  for  the  surface  area  per  unit  volume  of  the  atmosphere  in  the  plume. 

Atmospheric  measurements  of  the  alumina  particles  have  been  done  by  Strand,  et.  al.  (Ref. 
13),  and  by  Cofer,  et.  al.,  (Refs.  14, 15),  with  some  recent  laboratory  studies  by  Kim,  et.  al.  (Ref.  16). 
A  detailed  discussion  of  the  particle  size  distributions  has  been  recently  given  by  Beiting  (Ref.  17). 
The  work  by  Beiting  includes  the  possibility  of  a  variation  of  the  size  distribution  with  position  in  the 
plume.  In  our  calculations,  we  are  not  including  this  effect. 

The  measurements  by  Strand,  et.  al.,  were  done  at  20  km  altitude,  and  the  measurements  by 
Cofer,  et.  al.,  at  altitudes  up  to  7  km.  Gofer’s  work  clearly  shows  that  the  particles  are  spherical,  so 
the  surface  area  of  each  particle  is  simply: 


a  =  tiD^ 

where  D  is  the  diameter  of  the  particle.  For  a  given  size  distribution  dN/dD,  the  total  number,  N,  the 
total  area.  A,  and  the  total  volume,  V,  will  be  given  respectively  by: 

N  =  /  (dN/dD)  dD 
A  =  JjtD^  (dN/dD)  dD 
V  =  J(7t/6)D^  (dN/dD)  dD 


where  the  integrals  are  over  the  size  limits  of  the  distribution  function.  Examination  of  the  literature 
data  on  the  size  distribution  shows  that  a  single  function  for  dN/dD  is  not  adequate  to  represent  the 
data.  Instead  we  will  use  three  distributions  for  three  different  size  ranges  covered  by  the  data.  The 


A-1 


work  by  Strand,  et.  al.,  covers  sizes  from  0.02  to  3  microns,  and  requires  two  distribution  functions. 
The  work  by  Cofer,  et.  al.,  covers  roughly  1  to  10  microns,  and  requires  a  third  function.  We  will 
arbitrarily  match  the  Cofer  and  Strand  distributions  at  1.0  micron.  The  discussion  by  Beiting  (Ref. 
17)  gives  arguments  justifying  a  match  of  the  two  papers  at  this  point.  For  convenience,  we  will  put 
everything  in  terms  of  simple  exponential  distributions.  To  do  this  required  re-plotting  and  fitting  the 
raw  data  of  Strand,  et.  al.,  and  this  is  shown  in  Figure  14.  The  distribution  of  Cofer,  et.  al.,  is  used  as 
quoted  in  the  paper  (noting  a  typographical  error  in  the  paper.)  The  composite  fit  to  the  data  is  as 
follows: 

0.02  to  0.30  microns  diameter: 

dN/dD  =  4.1  X  10^ 

0.30  to  1.0  microns  diameter: 


dN/dD  =  18.4  e'^-^ 

1.0  to  10  microns  diameter: 
dN/dD  =  1  e 

The  composite  distribution  has  the  following  properties  when  integrated  over  the  respective  size 
intervals  (normalized  to  one  large  particle): 

Small  mode:  Nj  =  2.21  x  lO"*,  Dj  =  0.025  microns  (80%  of  surface  area) 

Medium  mode:  N2  =  3.64,  D2  =  0.552  microns  (3%  of  surface  area) 

Large  mode:  N3  =  1,  D3  =  2.39  microns  (17%  of  surface  area) 

The  literature  suggests  a  value  for  the  density  of  the  particles  of  2.0  gm/cin  (instead  of  4.0 
for  pure  alumina).  This  gives  for  a  final  value  of  area-to-meiss  ratio: 

4  2  2 

A/M  =  4.4  X  10  cm  /gm,  or  4.4  meters  /gm 


One  final  step  is  necessary.  SURFACE  CHEMKIN  is  set  up  to  work  with  a  monodisperse 
aerosol,  so  we  replaced  the  above  distribution  with  a  monodisperse  distribution  of  pure  alumina  (p  = 
4.0)  with  an  equivalent  surface/mass  ratio.  This  aerosol  has  a  uniform  diameter  of: 

D  =  0.34  microns. 
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